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ABSTRACT 

It is not yet clear what triggers the activity of active galactic nuclei (AGNs), 
but galaxy merging has been suspected to be one of the main mechanisms fuelling 
the activity. Using deep optical images taken at various ground-based telescopes, 
we investigate the fraction of galaxy mergers in 39 luminous AGNs (M^j < -22.6 
mag) at z < 0.3 (a median redshift of 0.155), of which the host galaxies are 
generally considered as early-type galaxies. Through visual inspection of the 
images, we hnd that 17 of 39 AGN host galaxies (43.6%) show the evidence for 
current or past mergers like tidal tails, shells, and disturbed morphology. In order 
to see if this fraction is abnormally high, we also examined the merging fraction of 
normal early-type galaxies in the Sloan Digital Sky Survey (SDSS) Strip 82 data 
(a median redshift of 0.04), of which the surface-brightness limit is comparable 
to our imaging data. To correct for the effects related to the redshift difference 
of the two samples, we performed an image simulation by putting a bright point 
source as an artihcial AGN in the images of SDSS early-type galaxies and placing 
them onto the redshifts of AGNs. The merging fraction in this realistic sample 
of simulated AGNs is only ~ 5 — 15% (1/4 to 1/8 of that of real AGNs). Our 
result strongly suggests that luminous AGN activity is associated with galaxy 
merging. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: general — 
galaxies: evolution — galaxies: interaction — galaxies: quasars 
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1. INTRODUCTION 


An active galactic nucleus (AGN) is an energetic core of a galaxy whose energy comes 
from accretion of surrounding matter by a supermassive black hole (SMBH) that resides 
at the central part of the galaxy. In order for an SMBH to shine as an AGN, it needs a 
supply of gas to fuel its activity. Two main mechanisms have been suggested to trigger 
AGN activity: an internal mechanism through a dynamical instability inside a galaxy and 
an external mechanism through galaxy-galaxy interaction or merging. However, it is not 
yet clear which one is the dominant mechanism, even after many observational studies have 
been carried out. 


The internal mechanism is such that a gas inflow to the central part occurs as a result 
of instability in the internal structure of a galaxy. For example, a galaxy bar can move gas 
from the outer regions of a galaxy into its center, and then the gas inflow can trigger the 
AGN phase 


( 

Lvnden-Bell 

1979: 

Sellwood 

1981; 

van Albac 

a & Roberts 

1981; 

Gombes & Gerin 

1985; 

Pfem 

lieer & Friedl 

1991 

; Heller & Shlosman 

1994 

Bournaud & Gombes 

2002; 

Athanassoula 


2003 


Sakamoto et ah 


19991). Studies of AGNs with moderate X-ray luminosity, especially 


the studies targeting distant AGNs, hnd that most of the AGN hosts are disk galaxies and do 


the AGN activitv and the SMBH growth ( 

Schawinski et ah 

2012. 2011: 

Gisternas et ah 

2011 

Kocevski et ah 2012; 

Simmons et ah 

2012). In the local 

universe, the fraction of Seyfert 


galaxies is found to be high in barred galaxies in the Sloan Digi tal Sky Survey (SDSS) data, 
suggesting a close connection between bar and AGN activ ity fjOh et al.l 120121) . but studies 
exist that point to an opposite conclusion flLee et al.l 120121) . 


In contrary, the external mechanism is represented by galaxy interaction or merg- 


Barnes & Hernauist 

1991; 

Mihos & Hernauistl 

1996 

Honkins et al.ll2005 

Snringel et al.ll2005) 


There are a number of observational results that support this idea. Studies of g alaxy pairs or 


galaxies in interaction hnd that the AGN fracti on increases in such systems. flEllison et ah 


201ll:ISilverman et al.lDOllt iBessiere et al.ll2012[) . Ultraluminous infrared gala xies (ULIRGs) 
are suggested to harbor AGNs, and many ULIRGs are merging systems (e.g., ISanders et ah 
1988al: I ll988bl) . Binar y SMBHs in sorne AGNs demon strate that two or more SMBHs can 


merge into one SMBH flKoss et al.ll201ll: lKomossall20031) . After all, many AGN host galaxies 
are found to be elliptical galaxies, which do not possess bars or disk instabilities. 


One promising way to investigate the AGN and merger connection is to study objects 
with merging features. When two galaxies with comparable mass merge, the merging pro- 
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duces an early-type galaxy flToomre &: Toomrel Il972l: iBarnes &: HernquistI Il992h. leaving a 
trace of the past merging activity in the form of ti dal tails, shells, and dust lanes flQninnlll 984J: 
Hernquist fc Spergellll992l: IFeldmann et alJl2008[l . In support of this theoretical expectation, 
very deep imaging of early-type galaxies hnd merging features in many cases (15% — 80% 


depending on the depth of the im age; Ivan DokkumI l2005l: iTal et al.l l2009l: iKavirail 12010 


Kim et ah 

2012: 

Sheen et ah 

2012) 


Feldmann et ah 


merging 

in so me models at the same timescale of roughly 1 Gyr after the merging flHopkins et ah 


2005h . Therefore, if major merging triggers AGN activity, we expect the following to be seen: 
(1) the AGN host will be an early-type galaxy or a bulge-dominated system, (2) merging 
features such as tidal tails and shells should be observable preferentially in AGNs than in 
nonactive early-type galaxies, and (3) AGNs should be found more in spheroidal postmergers 
(1-2 Gyr after the merging) than in merging systems at a very early stage. 

Interestingl y, previous sea r ches for merging features in AGN host galaxies have seen 
mixed results. iBennert et al.l (120081 ) found that four out of five low-redshift QSO hosts 
show obvious merging features from deep Hubble Space Telescope (HST) imaging data. A 
more extensive sample of 69 QSO hosts (0.142 < z < 0.198) was imaged by a ground-based 


telesc o pe, and it was fou nd that 60% of the sample shows signs of interaction fiLetawe et ah 


2010 |). 


Koss et ahl (120111) studied the fraction of mergers in the Swift BAT-selected AGNs and 


found that the fraction of BAT AGNs showing evidence of merger is much higher than the 
fraction in a sample of i nacti ve galaxies (18% versus 1%). For a special class of dust-reddened 
AGNs, lUrrutia et al.l (120081) also f ound a high fraction of AG N hosts with evidence for 
ongoing or recent interaction (85%). iRamos Almeida et al.l (120121) investigated the brightness 
of merging features in radio galaxies at 0.05 < z < 0.7 and found that the me rging features of 


radio-loud active galaxies are 2 mag brighter than those of quiescent galaxies. iGarpineti et ah 
(120121) found that about 70% of spheroidal postmergers —spheroidal galaxies with merging 
features —are Seyferts or LINERs. By comparing this fraction to ongoing merging systems, 
they suggested that the AGN phase becomes more d ominant only in the v ery hnal stage of 
the merging. On the other hand, an early study of iBahcall et al.l (jl997l) found that only 
15% of AGN hosts show a clear indication of interaction, although they also noted a large 
fraction AGN hosts with close, projected companions. However, the more recent studies of 
Karouzos et ah (2014a, 2014b) show that there is no enhancement of close neighbors around 
X-ray or radio-selected AGNs with respect to field galaxies with properties similar to the 
AGN hosts. The aforementioned studies supporting the internal AGN triggering mechanism 
all find that the fractions of both early-type host galaxies and merging systems are low in 
distant AGNs. 


These pieces of seemingly contradictory results can be resolved if AGN activities are 
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triggered by different mechanisms at different levels. There are theoretical studies that 
suggest s uch an idea, with the external mechanism preferentially triggering luminous AGN 
activity flHopkins fc HernquistI 120091: [Draper fc Ballantvnel l2012l ). From an observational 
point of view, the studies favoring the merging mechanism are mostly based on QSOs (bright 
AGNs) and early-type host galaxies. In comparison, the observational studies favoring the 
internal mech a nism are focused more on less-luminous AGNs or late-type host galaxies. 
Treister et ah! (12012! ) suggest that the fraction of merging systems in AGNs is a strong 
function of AGN luminosity, with merging activity being stronger in more luminous AGNs 
(AGNs with a large bolometric luminosity). Their sample, however, is rather heterogeneous, 
making it difficult to draw a hrm conclusion on this issue. Furthermore, a more recent 
study excludes a strong connection between AGNs and mergers over a wide range of X-ray 
luminosity (10^^ < L(X — ray)erg“^sec“^ < 10^^'®), based on their investigation of AGN 
host galaxies at 0 < 2 ; < 3 in the Ghandra Deep Field south, including a careful treatment 
of a control sample that is compared against the morpho logy of the AGN host g alaxies 
fjVillforth et al.ll2013|) . A similar conclusion was reached bv iGeoreakakis et ah! (12009!) too. 


Although researchers do not yet agree on the signihcance of merging for AGN triggering 
activity, it is known that host galaxies of lumin ous AGNs are predominantly early-type 
galaxies (iMcLure et al.l!l999l: [Dunlop et al.!!2003f) . i.e., products of major merging, so it is 


natural to consider merging as a key AGN triggering mechanism in luminous AGNs. The 
dehnition of “luminous” is a bit ambiguous, but the luminosity cut around ~ —23mag 
qualihes as “luminous” for the central point source. While previous studies have investigated 
the connection betwe en merging features and early-type hosts, the studies are rather limited 


galaxv tvpes and ar 

e not geared to detect the faintest features (e.g.. 

Georgakakis et ah 

2009; 

Gisternas et ah 

2011; 

Schawinski et ah 

2011 

Kocevski et ah 

2012 

Treister et ah 

20121. 


While deriving the fraction of merging signatures in AGNs can be useful, interpreting the 
result as a sign of merger-triggered AGN activity needs an additional step of comparison 
with an ordinary galaxy sample because many normal early-type galaxies are found to have 
merging features. However, such a comparison is a complicated process because the surface- 
brightness limits often differ between an AGN sample and a normal galaxy sample, and the 
bright point sources in AGNs can make it more challenging to detect merging features. 


In order to address the specihc question of whether luminous AGNs are triggered by 
merging activity, we study the fraction of luminous AGNs showing merging features and 
compare it to the fraction of normal galaxies with merging features. Because the light from 
luminous AGNs can contaminate the outer regions of their host galaxies, deep images of 
luminous AGN are obtained to reveal merging features. For a similar reason, we use a 
control sample of normal early-type galaxies selected from the SDSS Stripe 82 data, which 
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reaches a surface-brightness limit of Hr ~ 27 mag arcsec“^ over a 1" x 1" area at rms noise of 
Icr of the background regions of the imag^. In addition, merging features are more difficult 
to hnd in AGN hosts than in the galaxies in the control sample because of the redshift 
differences and the contamination from the central point sources. Therefore, we simulated 
images of luminous AGNs using galaxy images in the control sample so that the simulated 
AGNs contain a bright central point source and are placed at redshifts similar to those of our 
AGN sample. This enables us to do a realistic comparison between the control sample and the 
AGN sample. Finally, we compare the simulated AGNs versus the real AGNs, showing that 
merging features are much more easily found in host galaxies of luminous AGNs than in the 
simulated AGNs. Throughout the paper, we adopt a concordance cosmological parameters 
of Hq = 70 kms“^ Mpc, Qm = 0.3, and Ga = 0.7 (e.g., Im et ah 1997). The physical 
quantities taken from the literature were converted to match these cosmological parameters. 


2. AGN SAMPLE AND OBSERVATION 
2.1. Sample Selection of AGNs 

We constructed two samples of luminous type-1 AGNs at z < 0.3, one that we call as 
the “base sample” and another that we call as the “best sample,” which is a subsample of 
the base sample. The redshift cut of z < 0.3 is imposed because the features that we use to 
trace the merging activity are overmatched by the bright point source of the AGN or become 
too faint to detect if a galaxy is too far away. The definitions of “luminous” are different for 
the two samples, and they are explained in the following. 

The base sample comprises 39 luminous type-1 AGNs at z < 0.3. Here, the term 
“luminous” is defined to be Mr < —22.6 mag (Mr < —22.0 mag, or equivalently Mr < 
—22.6 mag assuming B-R = 0.6), and “type-l” AGNs are selected because they have an 
unambiguous spectroscopic signature (broad emission lines) for the existence of an AGN. 
The objects are chosen with this luminosity cut without separating the nuclear component 
from the host galaxy light, although we gave a preference to AGNs with bright nuclear 
components, such as Palomar Green (PG) quasars, during the sample selection process in 
order to select AGNs that are luminous in the nuclear component. Our analysis will focus on 
the base sample because the result from the best sample does not differ much from the result 
from the base sample, and the base sample provides better number statistics than does the 


^We will use this definitio n of SB l imit throughout the p aper. Even though less ri gorous than the definition 
of the SB limits in Li et al. ( 2011 1. Ho et al. ( 2011 1. and Bakos fc Truiillol ( 2012ll . this definition captures 
the depth of an image well in terms of how effectively we can identify faint, merging features. 
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best sample. The results from the best sample will be provided for the main results (Table 
5 and Figures 18 and 19) to demonstrate that the results from both samples are identical. 

We also constructed a subsample of AGNs (33) based on the nuclear magnitude cut, 
which we call as the “best sample.” We have good estimates of nuclear magnitudes from 
either HST or Canada France Hawaii Telescope (CFHT) images for a subsample of the base 
sample (59% or 23 out of 39). For the remaining 41% of AGNs, we estimated their nuclear 
magnitudes using a scaling relation between the black hole mass and the host galaxy bulge 
luminosity, as described in Section 1^751 For AGNs with good estimates of nuclear magnitudes 
(59% of the base sample), we hnd that most of them are dominated by a nuclear component 
with ruhost — ’^Nuc ~ 1 mag or tunuc ~ mtotai + 0.36 mag in median (see Section 13.2p . 
The absolute magnitude cut of Mr < —22.6 corresponds roughly to Mnuc(-R) < —22.24, 
assuming the above relation of uinuc ~ mtotai + 0.36 mag, and we adopt a slightly tighter 
nuclear magnitude cut of Mnuc(.R) < —22.44 to allow for a margin in the uncertainty in the 
nuclear magnitude estimates. The best sample will be used to support the validity of the 
results from the base sample. 


Note that AGN host galaxies are mostly early-type galaxies if their nuclear absolute 
magnitudes in the R band (Mnuc(-R)) are Mnuc(.R) < —22.44 mag. Figure [T] shows the 
histogram of nuclear magnitudes of AGN s divided into t wo types—-galax ies with B/T > 0.4 
which we consider as early-type galaxies film et al.ll2002l: iKim &: Imll2013l) . and galaxies with 
B/T < 0.4, which can be considered as late-type galaxies. Here, the data are taken from 
M. Kim et ah (2015, in preparation, hereafter, K15) where they decomposed the surface- 
brightness prohle of AGNs with the HST archival images into host galaxies and the central 
point source components. Figured] indicates that the host galaxies are dominated by early- 
type galaxies at Mnuc(-R) < —22.44 mag at the ~ 86.2% level. However, we did not impose 
the B/T > 0.4 cut in our AGN selection because the bulge absolute magnitude can be quite 
uncertain under the presence of a bright nuclear emission, with a typical error of ±0.2 - 0.4 
mag (Sanchez et ah 2004; Simmons et ah 2008; Kim et al. 2008b). The fraction of AGN 
hosts with B/T < 0.4 in our sample is small, anyway. 


The selection of AGNs was done in roughly two ways, depending on the facilities used, 
because the observations were performed at one of three observatories (Maidanak, McDonald, 
and Las Gampanas) or the data were taken from the archive (for two AGNs in our study). 


For observations at the Maidanak and McDonald observatories or for selecting objects 
from the archive, AGNs were selected from the Quasar s and Active Galactic Nuclei catalog of 


Veron-Gettv fc VeronI fj2010l ). Note that the catalog of IVeron-Gettv fc VeronI (1201011 provides 


the absolute B-band magnitude only. In order to convert the B-band ab solute magnitude to 


the R-band magnitude, we assumed a quasar composite spectrum from IVanden Berk et al. 
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fl200lh which gives B — R ^ 0.63mag. This value is between 0.46 mag, calculated from 
the slope ~ —1.03 in 6005 ~ 7180A derived by Glikman et ah ( 2006), and 0.71 mag, 
calcul ated from the slope a^, ~ —1.58 in 4000 ~ 8555 A derived bv IVanden Berk et ah 


(120011 ). We find 1,103 type-1 AGNs that satisfy the above sample-selection criteria. The 
absolute magnitudes listed in the Veron-Getty & Veron catalog, however, may include some 
amount of the light coming from the host galaxies. Therefore, we opted to observe the 
PG quasars among the possible targets that are observable during the assigned observing 
runs at the Maidanak and McDonald observatories because it is known that most of the 
PG quasars have a strong nuclear component (e.g., Kim et al. 2008a). At the end, the 
objects were selected mostly from the PG quasar catalog, except for two AGNs for which 
the data were taken from the GFHT archive. The selected PG quasars were inspected on 
the Digitized Sky Survey images, and none of them were found to be clearly dominated by a 
host galaxy component. Then, objects were randomly picked from a prepared target list for 
the observation. Priority was given during the observations to AGNs with lower redshifts 
because morphological features become more difficult to identify for higher redshift objects 
due to the surface brightness dimming and the resolution effects. In total, 18 AGNs of the 
sample and 14 AGNs of the best sample come from this selection. 


Th e sample observed with the DuPont 2.5-m telescope was taken from the iKim et al. 
(l2008al) and K15 samples, which are made of luminous type-1 AGNs with HST images. Kim 


et al. (2008a) and K15 assembled an imaging data set from the HST archive of type-1 AGNs 
for which signal-to-noise ratio values of the chosen HST images are high enough to disentangle 
the host galaxy comp onent from the br ight nucleus (Kim et al. 2008b). Using the HST 
images and GALFIT, iKim et al.l (j2008al) and K15 performed a two-dimensional, multiple- 
component surface-brightness profile fitting of AGNs, separating the nuclear component from 
the host galaxy and decomposing the host galaxy surface-brightness profile into a bulge, 
a disk, and a bar if necessary. The surface-brightness profile fitting provides structural 
parameters for the host galaxies. Typical uncertainties of B/T range from 0.05 to 0.2 (see 
also Ho & Kim 2014). Excluding three AGNs with host galaxies having a clear signature of 
spiral arms, the number of AGNs that are observed by the DuPont 2.5-m are 21 for the base 
sample and 19 for the best sample. 

Table [T] shows the list of AGNs in the base sample. Figure [2] shows the redshift versus 
the absolute magnitude of AGNs in the base sample in both Mr and These AGNs 

span a redshift range of 0.04 < z < 0.3, with a median redshift oi z = 0.155. 


In summary, there are 39 and 33 type-1 AGNs at z < 0.3 in the base and the best 
samples, respectively, where the best sample is based on a more rigorous sample-selection 
criteria than is the base sample, but at the expense of number statistics. Our base sample 




















includes all 13 PG quasars with Mr < —22.6 mag, 2 ; < 0.15, and decl. > 0 degree and a 
total of 23 PG quasars. 


2.2. Observation and Data Reduction 


In order to reveal faint, extended merging features around AGNs, we obtained deep op¬ 
tical images using several facilities. We observed 10 AGNs using SNUGAM (Im et al.2010; 
FOV= 18M X 18M and pixel scale= 0"266) on the 1.5m telescope at the Maidanak ob¬ 
servatory in Uzbekistan from 2010 June to 2011 August, six AGNs using the Gamera for 
QUasars in EArly uNiverse (GQUEAN; Park et ah 2012; Kim et ah 2011; Lim et ah 2013; 
FOV= 4'.8 X 4'.8 and pixel scale= Cf.'281) on the 2.1m telescope at the McDonald observatory 
in Texas from 2011 February to April, and 21 AGNs using the SITe2K GGD camera of the 
2.5m DuPont telescope (FOV= 8'.85 x 8'.85 and pixel scale= 0'.'259) at the Las Gampanas 
observatory from 2008 September to 2009 March. Table [U gives the summary of the hlter, 
the exposure time, the surface-brightness limit, and the observatory used for the observation 
of each target. The seeing ranges from 0'.'8 to 2" with a typical value at around 1". The 
adopted hlters range from V through i, depending on the observatory: The V band for 
the SNUGAM data, the i band for the GQUEAN data, the R band for the DuPont data, 
and the i band or the r band for GFHT archive data. We av oided using bluer hlte rs (e.g., 
B) because tidal debris are suggested to be rather red (e.g., iFeldmann et al.ll2008l ). Also, 
the choice of the hlter was driven by the observational efficiency to obtain deep images in 
minimal exposures. For example, GQUEAN at the McDonald observatory is optimized for 
red wavelengths and the r band suffers from aberration near the GGD edge, therefore the 
i band was chosen. For the Maidanak observation, the balance between faint sky and the 
desire to avoid a fringing pattern forced us to choose the V band. Although the chosen hlter 
sets are nonuniform, this is not likely to cause serious systematic errors in our analysis. 


Since it is easier to identify merging features at the fainter surface-brightness limit (/x), 
we tried to reach the surface-brightness limit of ~ 27 mag arcsec“^, comparable to the depth 
of the control sample. However, the resulting surface-brightness limit is a bit shallower 
than that of the control sample if the surface brightness dimming is taken into account for 
objects at z > 0.1. To reach a limit of ~ 27 mag arcsec“^, we observed each target with a 
series of short exposure frames with exposure times of 60-180s, avoiding saturation of the 
central point source in the case of the Maidanak images, and a combination of short (a few 
hundred seconds) and long (1,000-2,000s) exposures for the DuPont images. These frames 
are combined to create deep images with total exposure times between 1 and 3 hr, after 
properly rescaling images with different zero points and exposure times so that all of the 
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frames, before stacking, share the same zero point. In some cases, clear merging featnres 
became visible even at a short exposnre. In snch cases, the observation was stopped before 
reaching the desired snr face-bright ness limit (PG 0157-1-001 and PG 1613-1-658 at /r ~ 26.11 
and 26.18, respectively). The image of GTS J17.17 from the GFHT archive also shows a 
tidal tail at /i ~ 26.0 mag arcsec“^. 

Fignre[3]shows the redshift versns snrface-brightness limit of the base sample. Three dif¬ 
ferent symbols represent the three different filters (G, R or r, and i) nsed for the observation. 
Note that the snrface-brightness limit presented here is the “observed” snrface-brightness 
limit. The snrface brightness dimming canses the rest-frame snrface-brightness limit of each 
object to be brighter by ~ (1 -|- z)'^. This introdnces a bias in detecting merging featnres for 
an AGN at a higher redshift: onr AGN images are shallower than the control sample images 
for objects at a higher redshift. To indicate this effect, we also plotted with the dotted line 
the “effective” snrface-brightness limit of the control sample. This is a snrface-brightness 
limit of an image when galaxies in the image are assnmed to be at a certain redshift. For 
this line, we adopted the snrface-brightness limit of the SDSS Stripe 82 images as 27 mag 
arcsec”^ at their median redshift of 0.04. When these images are assnmed to be at a given 
redshift, the snrface-brightness limit at that redshift should be 27 -|- 10 log[(l -|- zj/l.Od] mag 
arcsec”^. As expected, we hnd that most of the AGN images are shallower than the dotted 
line or comparable to the surface-brightness depth within a few tenths of magnitude. We 
will make a correction for this bias through simulations of AGNs (Section 3.2). 


Basic reduction of the images, such as the bias and the dark subtraction, and the 
flat-fielding was done using the standard IRAF tasks. Flat-held images were taken during 
twilight, and these were used for hat holding the images. The dark correction was necessary 
only for the GQU EAN data where the dark current level is about 0.23 electrons/s/pixel 
(jPark et ahl 120121 1 and not negligible in the data with a long exposure time. The reduced 
frames were registered to the position of the hrst frame and stacked using the imcombine 
task of IRAF. Gosmic rays were rejected during the stacking process, and the bad pixels and 
columns were hxed by the hxpix task of IRAF before the registration of the frames. The 
photometry calibration was done by using standard stars observed for the SNUGAM and 
the DuPont data or by using the photometry information of stars from SDSS in the vicinity 
of the targets for the GQUEAN data. 


2.3. Black Hole Masses and Bulge, Nuclear, and Host Luminosities 

Because the merging fraction may depend on host-galaxy properties such as black mass 
or bulge luminosity, we provide these quantities for our AGN sample in Table [H 
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The black hole masses, Mbh in Table [T] are taken from iKim et al.l fl2008al) and 


Vestergaard &: PetersonI (120061) . where Mbh values are estimated by the reverberation 


mapping method for nine AGNs (marked as R in Table [T]) and by single-epoch methods for 
28 AGNs. Here, the single-epoch method refers to the methods that estimate Mbh values 
using spectra taken at a single epoch. In such methods, the Mbh values are derived from 
a combination of the 5100 A continuum luminosity and the line widths of broad emission 
lines (e.g., H/3 or Ha; Vestergaard & Peterson 2006). The Mbh for fhe remaining two AGNs 
are derived from the bulge magnitude of the host, which will be explained in a later part 
of this section. We caution readers that the absolute values of Mbh can be off by a certain 
amount. The Mbh estimates based on the virial method can be uncertain because they 
depend on the scaling factor (/), which is sensitive to the geometry of the broad line region. 
Recently, Kormendy & Ho (2013) showed that the correlation between the BH mass and 
the host galaxy mass in inactive galaxies is tight only for classical bulges and ellipticals. By 
comparing reverberation-mapped AGN to normal galaxies in the Mbh — o'* plane only for 
classical bulges and ellipticals. Ho & Kim (2014) recalibrated the scaling factor, which is 
approximately two times larger than the value used in this paper. Given the fact that our 
sample is hosted by early-type galaxies, it is very likely that bulges are either classical or 
ellipticals. Thus, it is possible that the BH mass is underestimated by a factor of two. 

The absolute magnitude of the bulge, Mbui, and the central point source, Mnuc, are 
estimated in two ways. In the hrst method, Mbui values are directly estimated from a 
two-dimensional surface-brightness htting that takes into account the nuclea r, bulge , and 
disk components. For AGNs with HST images, this is done in iKim et al.l (j2008al) and 
K15, and we take the results from their work. Note that the absolute magnitudes of the 
nuclear component and the bulge component are derived with K corrections that assume a 
power-law spectral energy distribution for the nuclear component (F)^ ~ z/" with a = —0.44 
blueward of 5, 000 A, and —2.45 otherwise; Vanden Berk et al. 2001) and an empirical 
template of elliptical galaxies for the bulge component. These a bsolute magnitudes include 
the correction for the Galactic extinction fjSchlegel et al.lll998l) . For the AGNs with the 
ground-based images only, we attempted this procedure for two AGNs with the archival 

resolved in these imag es. We htted the surface- 


GFHT images because host galaxies are wef 
brightness prohle of AGNs using GALFIT fjPeng et al.l 120021. l2010l) . assuming a model with 


three components, a point-spread function (PSF) for the central source and bulge and disk 
components for the host galaxy. The PSF was derived from a bright star in the vicinity of 
the AGN (< 6'). The disk component is r epresented b y an exponential prohle, and the bulge 
component is htted with a Sersic prohle (jSersid 119681) as given below: 



























11 


I{r) = Ieexp[-bn{r/reY^'^ - 1], 


( 1 ) 


where is the effective radius, Jg is the intensity at rg, and n is the Sersic index. During 
the fitting, the centers of each component were varied, and the background was fixed. We 
also fitted the object in the vicinity of the target to improve the fitting result. When the 
scale length of one component is found to be too small (< 1 pixel), we reduced the number 
of the models (PSF+bulge or PSF+disk) and used the deVaucouleurs profile for the bulge 
component. 


Figure m shows the CFHT images of CTS J17.17 and PKS 0159-062 fitted by GALFIT. 
The tidal tail and the jet-like features can be seen on not only the residual images but also 
in the original images of the two AGNs. Table [2] shows t he resultant fitting parameters. 
The absolute magnitudes are derived in the sa me manner as iKim et al.l fl2008al) . which uses 
an empirical template of an E/SO galaxy from iGalzetti et ahl (1199411 when deriving the K 
corrections. 


In the second method, Mr„i is estimated from Mbh using the correlation between Mbui 
and Mbh of active galaxies (IKim et al.ll2008al) as 


log(MBH/M0) — a + (3 Mbui{R), 


( 2 ) 


where a = —2.74 and (3 = —0.5. Then is obtained by subtracting this Mbui from 
the total luminosity of the AGN. This procedure assumes that Mbui ~ Muost and thus may 
overestimate the nuclear luminosity if there is a significant disk component in the host and 
underestimate the nuclear luminosity if the object is an outlier of Eq. (2). We used this 
method for six AGNs for whi ch neither HST nor GFHT irn ages were available. These are 
AGNs with Mbh values from Vestergaard fc Peterson (j2006 l. We used the relation of Kim 
et al. (2008) instead of the Mbh^Mbui correlation of inactive galaxies because host galaxies 
of AGNs often show signs of a young stellar population (e.g., Kauffmann et al. 2003; Ho 
& Kim 2014), leading to a somewhat smaller stellar mass-to-light ratio in comparison to 
inactive galaxies. 


Figure [5] shows the redshift versus the BH mass, the Eddington ratio, the bolometric 
luminosity, and the host bulge magnitude of our sample. Our sample spans parameter spaces 
of IO^-^Mq ~ IO^Mq, -2 < log(Lboi/-hEdd) < 0, 43 < log(Lboi/-h©) < 45.5, and -24 mag 
< Mbui( 7?) < —20.5mag. 
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3. CONTROL SAMPLE AND SIMULATION OF AGNs 


3.1. Control Sample 


We look into the merging fraction of AGNs through the presence of a merging feature 
in AGN host galaxies. To compare it with the merging fraction of inactive galaxies, we 
construct a control sample. The control sample is restricted to only early-type galaxies 
because the majority of host galaxies in our luminous AGN sample are early-type galaxies. 
For the control sample, we chose the SP SS Stripe 82 early-type galaxies with a surface- 
brightness limit of Pr = 27 mag arcsec“^ flKayiraill2010h. w hich is comparable to or deeper 
than the image depths of our AGN images (IKavirail l2010l) when taking into account the 
surface-brightness dimming effect in the images of the AGN sample. There are 317 early- 
type galaxies in this sample, which is limited to be at z < 0.05 and Mj. < —20.5mag. Table 
E] shows the basic properties of the control sample, such as the morphological type, B/T, 
and Mbh- 


These ear ly-type galaxies have been classihed through careful visual inspection by 
Kavirail 020 lOh . The classihcations include “relaxed,” “tidal feature,” “dust feature,” “tidal 
and dust feature,” “interacting,” and “spheroidal galaxy with faint disk features.” Among 
these classifications, we adopt “tidal feature,” “tidal & dust feature,” and “interacting types” 
as signs for a merging feature. There are 57 early-type galaxies identified as objects with 
merging features. A dust feature may be a sign of past merging, but we do not include it 
as a “merging feature” in our analysis because dust lanes do not extend to the outer region 
of galaxies and will be extremely difficult to detect in AGNs because a bright point source 
outshines the dust feature. Note that we simulate AGNs at different redshifts using the 
control sample in Section 3.2, in order to account for biases arising from a redshift difference 
between the control and the AGN samples, such as the surface brightness dimming and the 
change in the sizes of objects. A potential dependence of the merging fraction as a function 
of luminosity will be removed by analyzing the AGN sample at different host galaxy Mbh 
and luminosity bins and taking into account the difference in the Mbh distribution (Sections 
3.2 and 4 and Figure 9). 


The median redshift of the Stripe 82 early-type galaxies {zmed = 0.04) is lower than 
that of the AGN sample {zmed = 0.15). An evolution in the merging fraction could affect the 
analysis of the merging fraction difference between the control and the AGN samples, but 
such an effect is negligible (< 20%), as described below. Earlier works on the study of the 
merger fraction were controversial, with some studies showing a strong evolution in redshift 
that goes with ~ (1 -|- zY (e.g., Zepf & Koo 1989; Garlberg et ah 1994; Yee & Ellingson 
1995), and with others indicating a weak evolution (~ (1 -|- e.g., Neuschaefer et al. 
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1997). However, more recent results are converging to a rather weak redshift evolution in 
the merging fraction of (1 + a.t z < 1 (Lin et ah 2004; Lotz et al. 2008a; Bertone et ah 
2009; Bundy et al. 2009; Jogee et al. 2009; Lopez-Sanjuan et al. 2009; Man et al. 2012; Xu 
et al. 2012), which is especially true for massive or luminous galaxies that are comparable in 
the mass/luminosity range of the host galaxies studied here (Hsieh et al. 2008; Robaina et 
al. 2010; ). Even if we adopt the most extreme merging fraction evolution of (1 + z^ from 
the recent results, the expected increase in the merging fraction from ^ = 0.04 to z = 0.15 is 
only at the level of 20%, which is much smaller than the discrepancy in the merging fraction 
that we are exploring in this work. 


The merging features are further categorized by us into tidal tail (T), shell (S), or in¬ 
teraction (I), as shown in Figure O Our definition of “tidal tail” includes tidal tails and 
disturbed features like loops. The “shell” types are distinguished from “T” types as objects 
with distorted surface brightness and shell structures, but lacking features like loops. The 
“interaction” types are objects with their surface-brightness profile extending to a close com¬ 
panion. If there is no merging feature (i.e., “relaxed type” in Kaviraj 2010), we categorize 
such a case as no feature (N). Although there are quantitative and automatic classifica¬ 
tion techniques, we adopted the classification thr ough visual inspe c tion because the auto¬ 


matic classihcation may miss some faint features flLotz et al.ll2008bl: iKartaltepe et al.l 12010 
Adams et al. 2012). 


During the visual reinspection of the S PSS Str i pe 82 images, we found six galaxies 
classihed as objects with a merging feature in iKavirail (120lOl) . but we felt that it is not clear 
whether they really have merging features or not (Figure [7]). We treat them in our analysis 
as early-type galaxies with a merging feature, although we are unlikely to classify AGNs 
with this kind of object as AGNs with a merging feature. 


The bulge magnitudes of the control sample early-type galaxies are taken from the 
catalog of ISimard et al.l (1201 ll) . who performed a bulge-|-disk decomposition of galaxies in 
the SDSS DR7. We adopted the htted parameters from a model with free Sersic index n 
for the bulge compone nt and n = 1 for the disk component. This catalog covers about 84% 
of the control sample. ISimard et al.l (1201 ll) did not perform the bulge-fdisk htting for the 


remaining 16% of galaxies because they are brighter than mpetro(^) < 14 mag (too large to 
fit with within their automated pipeline). For these bright early-type galaxies, we assume 
that they are bulge-dominated galaxies. Note that the SDSS r-band magnitude is converted 
in to the R-band magn itude (r — i? = 0.151 mag), assuming an elliptical galaxy spectrum 
of iGalzetti et al.l fll994l ). The black hole mas ses were estimated using Eq. (2), but with the 
a = —2.6 and (3 = —0.5 of inactive galaxies flBettoni et al.l 120031) . 


Figure |8] shows the distribution of the bulge absolute magnitudes in R and the black 
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hole masses of the control sample as a function of 2 ;. Overplotted in the figures is our AGN 
sample. Figure |9] compares the Mbh distribution of the control sample and that of the 
AGN sample. There is a difference in the overall parameter space covered by the AGN and 
the control samples, where the Mbh values of the AGN base sample tend to be brighter 
and more massive than for the control sample. This difference can introduce a bias in the 
comparison of merging fractions if the merging fraction changes as a function of host galaxy 
brightness or black hole mass. Figure fTOl shows the merging fraction as a function of Mbh and 
Mbui of early-type galaxies in the control sample. The error bars indicate 68.3% confidence 
levels, which are computed using the method described in Gameron (2011)@. We hnd that 
the merging fraction increases with Mbh or the bulge luminosity. Therefore, we will take 
into account the difference in the parameter space (Mbh or Mbui) between the AGN and the 
control samples in our analysis, and we will use Mbh as the reference parameter to adjust for 
the difference in the covered parameter space. We adopt Mbh instead of Mbui(.R), because 
it is a better defined quantity for the AGN sample, and the Mbh versus Mbui relation is 
established more securely for passive galaxies (the control sample) than for AGNs. 


3.2. Simulation of AGNs Using a Control Sample 


It is well known that most (80 %) of early-type gal axies show merging features in very 


2P 

ml 


deep images (/i ~ 28 mag arcsec”^; Ivan Dokkumll2005l) . so a careful analysis is needed to 
analyze AGNs that are associated with merging . The fraction of early-type galaxies with a 
merging feature is 17% for the control sample fjKaviraill2010l) . but there are many factors 
that can influence the detection rate of a merging feature in the AGN sample. 


First, the depth of the image is important because many early-type galaxies with merg¬ 
ing features do not reveal merging features in a shallower image. We tried to match the 
depth of our AGN images and the control sample image, but not perfectly. Most notably, 
the rest-frame (effective) surface-bright limit of the AGN sample is a bit shallower by 0.5 - 1 
mag arcsec”^ at z = 0.2-0.3 than that of the control sample (^^med ~ 0.04) due to the surface 
brightness dimming. 

Second, the AGNs are located farther away than the control sample, so their apparent 
sizes are smaller. This causes a loss in the resolution element, which can wash away some of 
the merging features. This resolution effect also needs to be treated properly. 


Third, AGNs have bright nuclear components. For AGNs at a very low redshift, the 


^We will use this method to compute errors of fractional quantities throughout the paper. 
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PSF widths are small compared to the extent of the host galaxy, so the the bright nuclear 
component does not hamper the detection of a merging feature. However, if they are located 
at a higher redshift, the width of the PSF becomes comparable to or larger than the host 
galaxy size, and when combined with an extended wing of the PSF in many ground-based 
images, low surface-brightness features in the host galaxy becomes difficult to detect. 

In order to properly handle these biases, we perform a simulation where we take the 
control sample early-type galaxy images, plant a nuclear component, and place them at 
various redshifts. The simulation is done in the following way. 


In the hrst step, early-type galaxy images are block averaged by a factor of 2-5, which 
corresponds to the change in the angular diameter distance Da{z) from = 0.041 (the 
median redshift of the control sample) to the redshifts of 0.086, 0.138, 0.2, and 0.275, respec¬ 
tively, matching the redshifts of the AGN sample. The pixel scales of the block-averaged 
images are all considered as 073961/pixel, as in the original Stripe 82 image, to properly 
reflect the redshift changes. 


Because the simulated images are block averaged, FWHM of the PSF given in pixel 
units has shrunk. To mimic the PSF resolution (in arcseconds) of the AGN images, we 
convolved the block-averaged images using a Gaussian function with a ac value as 


<^0 “ ‘^out 


<^S82x2 

/ ^ 


where (Ts 82 is FWHM/2.35 of the image before block averaging, / is the block-averaging 
factor, and (Tout is the desired resolution of the output image after the convolution. We 
adopt two values of (Tout so that the seeing of the resultant image becomes either 1" or 2" 
in FWHM to allow us to investigate the effect of seeing in the identihcation of a merging 
feature. 


The pixel values are adjusted so that the total flux of the simulated galaxy is equal 
to the flux of the galaxy at the desired redshift. During the block-averaging process, the 
background noise is smoothed out. To make the surface-brightness limit of the resultant 
image to be ~ 27 mag arcsec“^, additional noises are implemented in the simulated images. 

Finally, we plant a point source at the center of the simulated early-type galaxy images 
in order to mimic an AGN and investigate the effect that a bright AGN component has 
on the merging fraction analysis. We hnd that the median of the differences between the 
host magnitude, rrihost, and the nuclear magnitude, rriNuc, is (nihost — 'nrNuc)med 1-0 mag, 
and the upper and the lower quartiles of the rrihost — values from the mode of the 

^^host — ■'^Nuc distribution are ~ 1.5, 0 mag. To simulate the nuclear point source, we chose 
a bright star in the SDSS Stripe 82 image and rescaled its pixel values so that the PSF 
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has three representative values of mgai — ruNuc that are 0, 1.0, and 1.5 mags brighter than 
the host magnitude. These PSFs are put at the center of the simulated early-type galaxy 
images. 

In summary, we made simulated AGN images for each of 57 SDSS Stripe 82 early- 
type galaxies with merging features, in a parameter that covers z = (0.041, 0.086, 0.138, 0.2, 
0.275), ruNuc = (none, rugai, nigai —1, rUgai —1.5), and seeing= (1", 2"), leading to the creation 
of 40 sets of simulated AGN images per a SDSS Stripe 82 early-type galaxy. Figures [TTl 
[Eiiinio show simulated AGN images based on three SDSS Stripe 82 early-type galaxies. 
Figures [TTl [13] are for the seeing of 1", and Figures [121 [HI are for the seeing of 2". Redshifts 
are 0.041, 0.086, 0.138, 0.2, and 0.275 from top to bottom, and the magnitudes of the nuclear 
point sources are none, mgai, rrigai-l.O, and mgai-1.5 mags from left to right. 

Figure [HI shows quantitatively how the merging features like tidal tails and shells become 
difficult to identify when the redshifts are changed and the nuclear sources are added. The 
top panel is for the case that does not include the nuclear point source, and it shows how 
the redshift effect alone affects the identification of a merging feature. The bottom panel in 
the figure assumes the addition of a nuclear source with mgai — 1 mag and a seeing of 1" and 
shows how the nuclear source makes it more difficult to see the merging features. We find 
that there is a steady decrease in the identihcation of a merging feature as we move objects 
to higher redshifts. With respect to z = 0.041, only 50% of the features can be identified at 
2 ; = 0.14, and this value goes down to 5 - 10% at z = 0.2. The shell features are more difficult 
to identify than the tidal features, and all of the shell features are virtually unidentihable at 
2 ; > 0.2. If we include a bright point source, the identification rate of merging features drops 
by 10 - 20% (or by a factor of two at an intermediate redshift). Our simulation demonstrates 
that the merging features can disappear easily at a higher redshift and with an addition 
of a bright nuclear component. Hence, one must take these effects into consideration when 
comparing the merging fraction of AGNs and with of the control sample galaxies. 


4. RESULTS 

4.1. Merging Features in the AGN Sample 

Figure [16] shows the images of our AGN sample in the order as listed in Table [T] Using 
Figure [T6l we identified the merging features through visual inspection. The merging features 
are categorized into tidal tail (T), shell (S), or interaction (I), and no feature (N), as was 
done for the control sample. The classification of the merging features is indicated in Table 
[Hand also in Figure [16] 


17 


Our visual inspection reveals that 17 of 39 AGNs (43.6%) show a merging feature. 
Among these AGNs with a merging feature, only one AGN shows a shell structure, and 
another shows a sign of interaction. Therefore, the majority of AGNs showing merging 
features are categorized as T. 


4.2. Merging Fraction: Comparison with Simnlated AGNs 


In the previous section, we hnd that 17 of 39 AGNs 43.6%) show merging features, 
but only ~17% of the control sample shows evidence of galaxy mergers. At face value, the 
merging fraction of AGNs is ~2.6 times higher than the merging fraction of the control 
sample. However, the difference in the merging fraction should be more pronounced if we 
consider other effects that could make it difficult to identify merging features in the AGN 
sample. For example. Figures [11] - [H] indicate that it becomes very difficult to identify 
merging features when an object is at a higher redshift and has a brighter nuclear component. 
We will take these effects into account and derive a more realistic comparison of the merging 
fraction between the control and the AGN samples. On the other hand, there are several 
AGN images that are deeper than the control sample images (the points above the dotted 
line in Figure 3). For these cases, the noises are added to the AGN images to match the 
control sample image depth, and the merging feature is classihed on the noise-added images. 
We hnd that the merging features can be identihed in the noise added images, which is 
not surprising considering that the depths of these AGN images are only a few tenths of 
magnitude deeper. 


Because the merging fraction changes as a function of Mbh (or equivalently, Mbui( 7?)), 
as was shown for the control sample in Figure 10, it is important that we use the control 
galaxies that have the same properties as the galaxies that AGN hosts evolve into. One could 
use Mbui( 7?) to tag control galaxies equivalent to an AGN host, but this is not desirable 


Woo et al. 

2006, 

2008; 

Hiner et al. 

2012 ) 


Moreover, 

the measurement of Mbui{R) is challenging for AGNs with a luminous nuclear component 
(Sanchez et ah 2004; Simmons et al. 2008; Kim et ah 2008b). Using the Mbh of an AGN 
host to identify its counterpart in the control sample is advantageous in this regard because 
BHs seem to have already grown at z ~ 0.3 in active galaxies, and Mbh estimators are 
well established for both AGNs and quiescent galaxies. However, the use of Mbh is not 
completely free f rom difficulties: the intrinsic scatter can be as large as 0.4 dex for Mbh 
estimators (e.g., ( Peterson et al.ll2004l) j. and BHs could still grow in high Fboi/FEdd AGNs 
(Kim et al. 2008). For the following analysis, we consider Mbh as a more desirable quantity 
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to match AGNs with equivalent galaxies in the control sample, but we also present results 
that use Mbui(-R) as a galaxy identifier in order to check if the choice of the galaxy identifier 
can change the main conclusion. 

Figure [T7] shows the merging fraction of the simulated AGNs as a function of redshift. We 
divide the sample into two, one with log(MBH/Af©) > 8.39 and another with log(MBH/Af 0 ) < 
8.39, because the merging fraction changes as a function of Mbh as we have shown earlier 
fSection l3.ip . The panels for each Mbh are also divided into three, each corresponding to 
a different nuclear magnitude. The magnitudes of the PSF components for each panel are 
’^gai (top), mgai — 1.0 (middle), and mgai —1.5 (bottom). The black asterisks are the merging 
fraction of early-type galaxies without any nuclear component included, and the black open 
circles are for the cases with a nuclear component included. The seeing of the simulated 
images is assumed to be 1" here. Tables IU [5] summarize the merging fractions in the figure, 
as well as similar numbers for the case of seeing at ~ 2 ". 

Figure fTTI shows that the merging fraction of the simulated AGNs decreases as a function 
of redshift, even without adding a nuclear component. At z = 0.041, it is 17%, but it 
decreases to 10% at z = 0.138 (a drop of 40%) and to 2% at z = 0.275. The effect of 
the addition of the nuclear component becomes prominent when the nuclear component is 
’^Nuc = ’^gai — 1 mag or brighter and at z > 0.138. At ^ > 0 . 2 , the merging fraction drops by 
50% or more due to the addition of the nuclear component. A larger seeing makes it more 
difficult to identify merging features. At z < 0 . 2 , the effect of the seeing is not much, with a 
decrease in the merging fraction by 20-30%. But at z > 0.2, it could decrease the merging 
fraction by 50% or more. 

In Figure [T71 we plot the merging fraction of our AGN sample too (the red stars). 
The AGN sample is divided into two redshift bins at ^ = 0.155 to provide enough number 
statistics in each bin. The horizontal bars on each point show the redshift range covered by 
the AGNs in each bin. The locations of the red points show the median redshift in each bin. 
We find that about 40 - 50% of AGNs with Mbh > 10 ®-^® Mq have merging features regardless 
of redshift. On the other hand, the merging fraction of AGNs with Mbh < 10 ®'^® Mq varies 
from 60% to 20% going from low to high redshift. Note that the actual merging fraction 
of AGNs could be higher because the redshift shift and the bright nuclear component can 
make the merging feature difficult to identify, as we saw from the simulation in the previous 
section. In comparison to the simulated AGNs, the merging fraction of AGNs is much higher 
by a factor of 3-10, except for the redshift bin at z = 0.12 of AGNs with massive black holes. 

The result is combined into Figure (THl where we plot the merging fraction of AGNs 
at two different black hole mass bins, and compare it to the same quantity derived for 
the simulated AGNs (the seeing value of 1" is assumed). The red stars show the merging 
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fraction of AGNs in the base sample in two mass bins (7.3 < log(MBH/Af 0 ) < 8.39 and 
8.39 < log(MBH/Af 0 ) < 9.0). The red hlled stars show the merging fraction of AGNs in the 
best sample. The black hlled circles are the merging fraction of the simulated AGNs with 
’^Nuc = ’^gai — 1 mag, i.e., the median of the nuclear magnitudes of the AGN sample, and 
the black points are plotted at the median Mbh of each Mbh bin. When calculating the 
merging fraction of the simulated AGNs, we use Eq.([3]) in order to take into account the 
difference in the Mbh distribution between the control sample and the AGN sample (Figure 
9), the spread in the redshift distribution of AGNs (Figure 5): 



,Sz{j) -^control,merg (5MBH,i, SZj) 


(3) 


Merging fraction 


fcorr,SMBii(i) -^control 
k 


Here, the quantity Ncontmi,merg{SMBu,i,Szj) is the number of simulated AGNs that are 
identihed to have merging features with their Mbh and redshifts belonging to an Mbh bin of 
5MBH,i (e.g., 8.25 < log(MBH/M 0 ) < 8.5) and a redshift bin of 5zj, and A^"controi(<^MBH,i) is 
the total number of control galaxies in the mass bin of ^Mbh.i- The factors fcorr,SMBH{i) 
fcorr, 5 z{j) correct for the difference in the Mbh distribution between the AGN and the control 
samples and account for a range of redshifts of AGNs. They are dehned as fcorr,SMBH{i) — 
[A^controi(AMBH)/A^controi((5MBH,i)] X [A^agn(^M nH.i)/^AGN(AM bh)] and 

fcorrMj) = NAGT<i{Szj, AMbh)/Aagn(AMbh), where AMbh represents either one of the 
two broad mass bins of 7.3 < log(MBH/M 0 ) < 8.39 and 8.39 < log(MBH/M 0 ) < 9.0. 

Irrespective of whether we adopt the base or the best sample, we hnd the following 
(Table [HD . At the higher mass bin of 8.39 < log(MBH/M 0 ) < 9.0, the AGN merging fraction 
is about 48 ± 15%, and it is 12 ± 5% for the simulated AGNs, which is a difference of a factor 
of four. At the lower mass bin of 7.3 < log(MBH/M 0 ) < 8.39, the AGN fraction is about 
39%, and that of the simulated AGNs is only 5 ± 2 %, giving a factor of eight difference. 

Similarly, we plot the merging fraction as a function of the absolute magnitude of the 
bulge (Figure [T9D. The same method is adopted for computing the merging fraction as for 
the Mnn-based merging fractions but using Mbui(7?) as a quantity to compute the correction 
factor. We hnd exactly the same trend as in Figure [T8l although the difference in the merging 
fraction between the control and the AGN samples is reduced in at the lower luminosity bin 
(with respect to the lower Mbh bin). The reason for the discrepancy can be understood 
when we examine the difference in the overall distribution of AGNs in Mbh versus with 
respect to the control sample, as seen in Figure [HI In Mr, most AGNs occupy the luminous 
end of the control sample, while it is less so in Mbh- Luminous early types in the control 
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sample have a higher merging fraction than fainter early types, which leads to the higher 
merging fraction at the lower luminosity bin of the control sample. However, the amount of 
luminosity evolution is uncertain for early-type host galaxies of AGNs, and the amount of 
luminosity di mming can be w’ell above a few tenths of t he magnitude of normal early-type 


galaxies (e.g., Ilm et al.lll996l. 120021: iBernardi et ahll2003h from = 0.3 to 0.04. If we make 


such a correction for the luminosity dimming in the AGN host galaxy, the Mr distribution 
would resemble the Mbh distribution more closely, and the merging fraction of the control 
galaxies at the lower luminosity bin would decrease to the level of the value for the lower 
Mbh bin. 


We note that our classihcation using visual inspection is not likely to affect our result. 
To quantify possible errors in the visual classihcation, the visual classihcation was done by 
two authors independently for the simulated AGN sample. We hnd that the classihcation 
agrees fairly well, and the resultant error in the merging fraction is only about ±0.01 to 
±0.03. Overall, this is a negligible amount. 

We also note that a slight diherence in the redshift ranges of the control sample and 
the AGN sample does not bias our result in favor of higher merging fraction in the AGN 
sample. As we already argued in Section 3.1, the merging fraction evolves only mildly over 
the redshift range oi z = 0.15 and = 0.04. The redshift evolution trend can be reduced by 
limiting our AGN sample to those at z < 0.1. Doing so only makes our result stronger, with 
hve out of six AGNs at z < 0.1 showing merging features. 

When dehning our sample, we excluded AGN hosts with spiral arms when HST images 
are available. We could not do so for AGNs whose images come from the ground-based 
observation. This can lead to a slight underestimate of the AGN merging fraction because 
AGN hosts with spiral arms (i.e., a significant disk component) tend to have small Mbh 
values and therefore a smaller merging fraction than those with large Mbh values. 

Although we hnd that 43.6% of luminous AGNs in our sample show merging features, 
our simulation of AGNs using the control sample strongly suggests that this is more likely a 
lower limit of the merging fraction. The true merging fraction of the luminous AGNs can be 
estimated by using the simulation result. We have shown earlier that the merging fraction 
is reduced by a factor of 1.7-7.6 at z = 0.15 to z = 0.2 with respect to the merging fraction 
at z = 0.041. If we make this correction to the observed merging fraction of AGNs, then the 
merging fraction becomes about 75-100%. Overall, we conclude that the merging is closely 
related to luminous AGN activity and very much likely to be the triggering mechanism of 
luminous AGNs. 
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DISCUSSION 


Tidal tails can be more prominent in gas-rich “wet mergi ng” fiFeldmann et al.N2008l) 
and appear from an early stage of galaxy merging fjMihoslll995l) . As the distnrbed gaseons 
and stellar materials fall back to the merger prodnct, the tidal tail shrinks and the fal l- 
back material can form shell-like or loop-like strnctnres flMihoslll995l: iFeldmann et al .11201)811 . 
Therefore, we snbdivided the merging featnre into the tidal-tail type (early stage of merging), 
the shell type (late stage of merging), and the interaction type (very early stage of the 
merging), as an attempt to see which stage of the merging onr AGN hosts are fonnd. We 
hnd that only one of 17 AGNs with a merging featnre shows the shell strnctnre. On the 
other hand, the majority of AGNs with merging features (15 out of 17) are in the category of 
the tidal tail type, and the remaining one is the interaction type. In comparison, ~ 25% of 
galaxy mergers in the control sample show the shell structure. At a hrst glance, this result 
appears to be in perfect agreement with the merger-triggered AGN scenario. However, we 
have shown through our simulation that shell structures are difficult to identify when objects 
are farther away a nd have a brig h t cent ral point source. For example, three of our AGNs 
are also studied in IBennert et al.l (l2008l) . and they show that two of them (MG 1635-1-119 
and PKS 0736-1-01) have a shell-type feature in their HST images, but we do not see such 
a feature in our low-resolution ground-based image, clearly indicating that we are missing 
the shell-type merging features in some cases. Therefore, we conclude that the lack of tidal- 
feature-type AGNs in our sample is consistent with a merger-triggered AGN scenario, but 
a hrm conclusion on this issue should be obtained through future studies with deeper and 


higher resolution images of AGNs (e.g.. iBennert et al.l 120081) . 


We mentioned earlier that we could not identify merging features f or some contro l 
galaxies that have been marked as early types with merging features in iKavirail (120101) . 
Gonsequently, such objects would not be classihed to have a merging feature even if they 
exist in the AGN sample. A correction to this difference in the merging feature classihcation 
between classihers (Kaviraj vs. us) could only more strongly in favor of our result. 

Because some of our AGNs overlap with those studied in previous works, we can make 
a direct comparison of the classihcation between ours and the oth er works. We hnd t hat 1 8 
objects in our AGN s ample overlap with th e AGN hosts studied by Hutchings fc Neh_(|l992 '). 
Bahcall et al. ( 1997 ). McLure et ah ( 1999 ). Dunlop et ah ( 2003 ). Bennert et ah ( 2008 ). and 
Let awe et ahl (120101) . These studies are based on either HST images or ground-based images. 
We hnd that three AGNs in the overlap sample are classihed to have merging features in 
our study (PKS 1020-103, PKS 2355-082, and PG 1302-102), but they are mentioned only 
as objects with companions in the other studies, demonstrating that our deep images newly 
uncovered faint merging features in some AGNs. On the other hand, we miss two AGNs with 
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a shell-type merging feature as mentioned above, showing that our images are not sensitive 
to picking up smooth, symmetric, low surface feat ures with detai l ed str uctures such as shells. 
One object, PG 0052-1-251, was classihed as Sb in lBahcall et ahl (119971) . but we classify it as 
a tidal-tail type. The s piral structure is mentioned to be terminated at a companion galaxy 
Dunlop et ahl (l2003l) . so the spiral structure may be related to a merging activity. The 


m 


bulge-fdisk decomposition of the host surface-brightness profile indicates that B/T = 0.67 
(K15), suggesting that the host galaxy is an early-type galaxy, lending a farther support 
that the spiral structure is related to the past merging activity. Our classification agrees 
with that of previous works for the remaining 12 AGNs in the overlap sample. 

Finally, we comment on the contradictory observational results regarding the AGN trig¬ 
gering mechanism. Our results that merging features are much more frequently found in lumi¬ 
nous AGNs than in quiescent galaxies, and that the AGN merging fraction could be as high as 
75%-100%, point toward a strong connection between merging and luminous AGN activity. 
Our AGN sample spans a range of lO^^er^fs”^ < Lboi < with a median bolomet- 

ric luminosity of Lboi ~ and the AGN samples of studies that have not found 

merging evidence in AGNs extend to a much fainter luminosity range (Tboi ~ lO'^^er^fs"^). 
Therefore, the difference in the explored luminos ity range could be the main source of the dis¬ 
crepan cy, as suggested in previous works such as Treister et ah fj2012 h Hopkins fc Hernquist 


(120091) . and (Draper fc Ballantvnd (120121) that major merging is only dominant in luminous 
AGNs. This leads us to a question whether fainter AGNs are really triggered by an in¬ 
ternal mechanism. After all, SMBHs are predominantly found in bulges of spiral galaxies 


20061 : iHopkins et ahl l2010l : IScannapieco et ah 


hint AGNs. 


that are the main hosts of 
produced through merging (lAeuerri et al 


however, bulges of late-type galaxies can also be 


2001 


2010 


Kannaupan et al.ll2004l: lEliche-Moral et al. 


Oser et al.ll2012h . Then, if the produc¬ 


tion of the bulges is closely connected to AGN activity, much like in the case of major 
merging-triggered AGNs, merging features from merging events (possibly minor merging) 
could be detectable in less luminous AGNs too. The tidal tails from faint satellite galaxies 
are going to be much less prominent than the tidal tails associated with major mergers in 
luminous AGNs, simply if the brightness of tidal tails scales with host galaxy luminosity. 
Thus, images that go deeper than current survey limits may reveal traces of the past minor 
merging activities in less luminous AGNs. 


For example, studies based on HST survey images have shown little evidence for merging 
activities in less luminous AGNs. Such imaging data, although very deep in detecting point 
sources, suffer in depth when the surface-brightness limit is considered. This is largely 
due to the surface brightness dimming effect that scales as ~ (1 -|- z)^. In the case of the 
GANDELS survey (Grogin et al. 2011; Koekemoer et al. 2011) and the images that were 
used by Schawinski et al. (2012) to draw a conclusion that only 4% of AGNs at z ~ 2 
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went through major merging, the SB limit is only 25.4 AB mag arcsec“^ at F775W (similar 
to R band) when the SB dimming correction is made. Future studies with deeper imaging 
data will certainly tell us whether the AGN triggering mechanism is closely related to the 
formation mechanism of bulges. If the classical bulges are formed through merging but the 
pseudobulges are formed through an internal mechanism (e.g., Kormendy & Ho 2013), such 
a difference will be reflected in the existence or nonexistence of merging features around 
luminous AGNs and in the photometric properties of bulges. 


6. SUMMARY 

To investigate whether galaxy merging can trigger AGN activity or not, we examined 
deep images (surface-brightness limit of /i of ~ 27 mag arcsec”^) of 39 luminous AGNs 
at Mr < —22.6, 2 : < 0.3, and 7.3 < log(MBH/Af 0 ) < 9.0. We hnd that 17 of 39 AGNs 
(~ 43.6%) show evidence of galaxy mergers such as tidal tail, shell, and interaction through 
careful visual inspection. Gompared to the control sample of early-type galaxies taken from 
the SDSS Stripe 82 data at a similar surface-brightness limit, the fraction of AGNs showing 
the merging feature is about 2.6 times higher. This difference becomes more signihcant when 
we simulate AGNs using the SDSS Stripe 82 early-type galaxy images because our AGNs 
are located farther away than the SDSS Stripe 82 galaxies (surface brightness dimming), 
and also bright nuclear components of AGNs can make it more difficult to detect merging 
features. Overall, we hnd that the merging fraction of the simulated AGNs is only 5-15%, 
much smaller than our result by factors of 4-8 depending on Mbh values. Our result strongly 
suggests that merging plays an important role in triggering AGN activity of luminous AGNs. 
Future studies with deep imaging observation of fainter AGNs are needed to understand how 
AGN activities are triggered in fainter AGNs. 
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Table 1. AGN Sample 


Name 
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-22.34 

-22.34 

- 

V 

10800 

27.15 

1.86 

Ma 

N 

PG 1402+261 

0.17 

15.34 

-24.57 

7.5 

S(Kim) 

-23.23 

-21.17 

-22.51 

0.486 

V 

12960 

27.06 

1.47 

Ma 

N 

PG 1411+442 

0.09 

14.01 

-24.57 

8.4 

R(Kim) 

-23.6 

-21.54 

-22.3 

0.076 

i 

12600 

27.40 

2.54 

Me 

T 

PG 1416-129 

0.13 

16.1 

-22.97 

8.5 

S(Kim) 

-22.82 

-21.59 

-21.59 

0.044 

R 

2400 

27.37 

0.92 

D 

N 

PG 1440+356 

0.08 

14.58 

-23.37 

7.4 

S(Kim) 

-23.19 

-20.91 

-22.17 

0.608 

V 

7550 

27.52 

1.17 

Ma 

T 

PG 1519+226 

0.14 

16.5 

-22.97 

7.9 

S(V&P) 

-22.66 

-21.38 

-21.38 

- 

i 

9000 

27.57 

2.08 

Me 

N 

PG 1613+658 

0.13 

15.49 

-23.37 

8.2 

R(Kim) 

-24.32 

-23.08 

-24.08 

0.297 

i 

5400 

26.18 

2.02 

Me 

T 

PG 1617+175 

0.11 

15.39 

-23.37 

8.7 

R(Kim) 

-23.22 

-21.21 

-21.38 

0.029 

R 

3600 

27.41 

0.88 

D 

N 

PG 1626+554 

0.13 

15.68 

-23.67 

8.5 

S(V&P) 

-23.22 

-22.5 

-22.5 

- 

i 

6000 

27.32 

2.54 

Me 

N 

PG 2214+139 

0.07 

14.66 

-22.96 

8.6 

S(V&P) 

-21.56 

-22.61 

-22.61 

- 

V 

10800 

27.36 

1.1 

Ma 

S 

PKS 0159-062 

0.19 

16.6 

-23.17 

8.7 

- 

-20.94 

-22.87 

-22.9 

- 

r 

2180 

27.96 

0.73 

C 

T 

PKS 0736+01 

0.19 

16.47 

-23.47 

8.1 

S(Kim) 

-24.12 

-23.39 

-23.39 

0.381 

R 

3600 

27.52 

1.48 

D 

N 

PKS 1020-103 

0.2 

16.11 

-24.17 

8.7 

S(Kim) 

-23.38 

-23.02 

-23.02 

0.061 

R 

3600 

27.20 

1.19 

D 

I 

PKS 1217+02 

0.24 

15.97 

-24.77 

8.3 

S(Kim) 

-24.28 

-23.4 

-23.4 

0.278 

R 

3600 

27.02 

0.74 

D 

N 

PKS 2135-14 

0.2 

15.53 

-24.87 

8.9 

S(Kim) 

-23.98 

-23.12 

-23.12 

0.043 

R 

4440 

26.48 

1.57 

D 

N 

PKS 2349-01 

0.17 

16.59 

-23.27 

8.6 

S(Kim) 

-23.81 

-23.34 

-23.75 

0.12 

R 

3080 

27.25 

0.84 

D 

T 

PKS 2355-082 

0.21 

17.5 

-23.07 

8.6 

S(Kim) 

-23.02 

-23.32 

-23.5 

0.034 

R 

4680 

26.85 

0.94 

D 

T 

Note. — Column 

1: object name; 

column 2: 

redshift; 

column 3: 

apparent 

magnitude 

in the V band 

from IVeron- 

•Cettv Veronl J201GI+ 

column 4; 

absolute 

J?-band 

magnitude, 


converted from absolute 5-band m agnitude; col umn 5: the black hole m a ss; co lumn 6: references for Mbh and method for estimating R=reverberation mapping; S=single-epoch 

method; Kim ^Kim et al.l i2008al) ; K15; V&:P= jVestergaard PetersonI ll2006l~) : column 7: absolute i?-band magnitude of the AGN nuclear component; column 8: absolute K-band 
magnitude of the bulge of the AGN; column 9: absolute R-band magnitude of the AGN host galaxy (see section 12.31 for how we derived quantities in columns 5 and 7—9). Bulge 
magnitude s come from an SB fitting of the images from D or C in the column 15 and from the black hole mass for Ma or Me in the column 15; column 10: Eddington ratio from 
iKim et al.l ll2008al~) : column 11: filter; column 12: integrated exposure time; column 13: surface-brightness limit, averaged over a l" X 1 " area at Icr per pixel in the observed band; 


















column 14: seeing FWHM in arcseconds; column 15: telescope: D = 2.5m DuPont telescope in Las Campanas, C = CFHT, Ma=1.5m telescope in Maidanak observatory, Mc=2.1m telescope 
in McDonald observatory; column 16: merging type T=tidal tail; S = shell structures; I=interaction; N = no merging feature. 

“ Typical formal error is <0.1 dex, but the error arising from the intrinsic scatter of the M-qh estimators could be as much as 0.4 dex. 

^Typical error is ~ 0.1 mag. 

‘^Typical error is zbO.3 — ±0.4 mag. 

‘^Typical error is ±0.2. 



Table 2. Fitted Parameters 


Name 

^Nuc 

mBui 

Re 

n 

'^disk 

Rd 

Hd 


Filter 


(mag) 

(mag) 

(") 


(mag) 

["] 




(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

CTS J17.17 

16.74 ± 0.00 

16.24 ± 0.01 

2.2 ± 0.01 

1.27 ± 0.02 

16.4 ± 0.01 

6.23 ± 0.16 

1 

1.18 

i 

PKS 0159-062 

19.11 ± 0.01 

17.37 ± 0.00 

4.36 ± 0.03 

[4] 

- 

- 

- 

1.11 

r 


Note. — column 1: object name of CFHT data fitted through GALFIT; column 2: apparent nuclear magnitude in the observed 
filter; column 3: apparent bulge magnitude in the observed filter; column 4: effective radius of the bulge; column 5: Sersic index of 
the bulge; column 6: apparent disk magnitude in the observed filter; column 7: disk scale length; column 8: sersic index of the disk; 
column 9: reduced column 10: filter. 








Table 3. Control Sample of Early-type Galaxies 


Sequence 

(Number) 

( 1 ) 

SDSS ID 

( 2 ) 

R.A. 

(deg) 

(3) 

Decl. 

(deg) 

(4) 

z 

(5) 

R 

(mag) 

( 6 ) 

Error 

(mag) 

(7) 

BjT 

( 8 ) 

Error 

(9) 

Mr 

(mag) 

( 10 ) 

Error 

(mag) 

( 11 ) 

Mbh 

log(M/M 0 ) 

( 12 ) 

Type 

Kaviraj (2010) 
(13) 

Type 

This work 
(14) 

2 

587730845817504263 

321.915192 

-1.102175 

0.030 

14.739 

0.003 

0.150 

0.000 

-18.800 

0.020 

6.875 

1 

N 

10 

587730846891507736 

322.370728 

-0.219635 

0.050 

14.415 

0.002 

0.410 

0.010 

-21.320 

0.020 

8.135 

4 

T 

13 

587730846891507870 

322.451935 

-0.313562 

0.030 

14.518 

0.002 

0.550 

0.000 

-20.490 

0.010 

7.720 

3 

N 

24 

587730847963939362 

319.503601 

0.570868 

0.035 

15.062 

0.003 

0.940 

0.020 

-20.980 

0.020 

7.965 

1 

N 

49 

587730848501596412 

321.287964 

0.896778 

0.049 

15.267 

0.003 

0.490 

0.010 

-20.660 

0.010 

7.805 

1 

N 

50 

587730848501399575 

320.744537 

1.026719 

0.031 

14.173 

0.002 

0.500 

0.010 

-20.630 

0.010 

7.790 

1 

N 

76 

587731174382567703 

320.190460 

0.790999 

0.043 

15.405 

0.003 

0.990 

0.010 

-21.080 

0.010 

8.015 

1 

N 

89 

587731185116053719 

327.825348 

-0.957031 

0.027 

13.665 

0.002 

1.000 

-99.990 

-21.584 

0.002 

8.267 

2 

T 

112 

587731185125818492 

350.121124 

-1.002414 

0.031 

13.519 

0.002 

1.000 

-99.990 

-22.015 

0.002 

8.483 

1 

N 

128 

587731185650565143 

322.302155 

-0.621743 

0.030 

14.433 

0.002 

0.520 

0.000 

-20.360 

0.010 

7.655 

2 

T 

129 

587731185650565167 

322.305756 

-0.523884 

0.030 

14.723 

0.003 

0.530 

0.000 

-20.060 

0.010 

7.505 

1 

N 

135 

587731185651155402 

323.751648 

-0.511436 

0.030 

14.439 

0.002 

0.500 

0.000 

-20.270 

0.010 

7.610 

2 

T 

146 

587731185654366563 

331.085999 

-0.502819 

0.046 

15.855 

0.004 

0.340 

0.010 

-19.310 

0.040 

7.130 

1 

N 

192 

587730846893080588 

325.963776 

-0.321398 

0.027 

14.169 

0.002 

0.570 

0.010 

-20.400 

0.010 

7.675 

1 

N 

221 

587731186724110554 

321.939697 

0.326403 

0.031 

14.692 

0.003 

0.750 

0.010 

-20.540 

0.020 

7.745 

1 

N 

237 

587731186735186207 

347.261871 

0.266889 

0.033 

15.056 

0.003 

0.560 

0.010 

-19.680 

0.020 

7.315 

2 

T 

261 

587731187271663782 

346.281830 

0.826253 

0.042 

15.627 

0.003 

0.960 

0.060 

-20.980 

0.070 

7.965 

2 

T 

266 

587731187279069315 

3.300832 

0.742940 

0.039 

14.342 

0.002 

0.880 

0.010 

-21.640 

0.010 

8.295 

1 

N 

269 

587731187279593490 

4.393614 

0.743043 

0.044 

14.596 

0.002 

0.880 

0.010 

-21.790 

0.010 

8.370 

2 

T 

272 

587731187279462423 

4.088767 

0.788624 

0.044 

15.253 

0.003 

0.700 

0.010 

-20.790 

0.010 

7.870 

1 

N 


Note. — column 1: sequence number of Kaviraj (2010); column 2: SDSS object ID; column 3: R.A. (J2000); column 4: Decl. (J2000); column 5: redshift; column 6-7: apparent J?-band 
magnitude and its error; column 8-9: BjT and its error. The error of -99.990 indicates that the SB model does not have a disk component; column 10-11: absolute K-band magnitude of 
the bulge component and its error; column 12: black hole mass; column 13: morphological classification from Kaviraj (2010): l:relaxed, 2:tidal feature, 3:dust feature, 4:tidal and dust 
feature, 5:interacting, 7:spheroids with faint disk; column 14: morphological classification from this work: N:relaxed, T:tidal feature, S:shell-type. (This table is available in its entirety in 
a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 
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TTlgal 

rrigai -1.0 

mgal -1.5 



loglMBH/AT©) < 8.39 



log(MBH/Af0) > i 

.086 

0.138 

0.2 

Redshift 

0.275 

0.041 

0.086 

0.138 

1.116 

0.074 

0.0141 

0.011 

0.455 

0.394 

0.303 

1.106 

0.067 

0.014 

0.007 

0.424 

0.364 

0.242 

1.095 

0.042 

0.004 

0.004 

0.394 

0.303 

0.212 

1.085 

0.042 

0.004 

0.004 

0.364 

0.242 

0.212 



Note. — The simulated AGNs are divided into two Mbh bins, and the results are for the seeing of l" 
They are simulated to be at several different redshifts and to have various nuclear magnitudes. 


Table 5. Merging Fraction of Simulated AGNs II 


1 // 2 " 
Redshift 


"^Nuc 

0.041 

0.086 

0.138 

0.2 

0.275 

0.041 

0.086 

0.138 

0.2 

0.275 

None 

0.167 

0.142 

0.098 

0.022 

0.016 

0.164 

0.132 

0.069 

0.009 

0.003 

?7T.gal 

0.167 

0.132 

0.085 

0.022 

0.013 

0.155 

0.12 

0.057 

0.009 

0.003 

mgai -1.0 

0.151 

0.11 

0.06 

0.009 

0.006 

0.145 

0.085 

0.041 

0.003 

0.00 

mgai -1.5 

0.151 

0.11 

0.06 

0.009 

0.006 

0.145 

0.085 

0.038 

0.003 

0.00 


Note. — The results for two seeing conditions (l" and 2") are presented for simulated AGNs placed at 
several different redshifts and with different nuclear magnitudes. 


Table 6. Merging Fraction of AGNs versus Simulated AGNs 




AGN 

Simulated AGN 

log(MBH/Af0) 

Base 

Best 

mpsF = rngai — 1 



(Mnuc(+ < -22.44) 

V seeing 

7.3 < log(MBH/MQ) < 8.39 

0.389t°;J2i 

n TTT +0-135 

U.OOO _o.096 

n 051 +0.017 

U.UOi _o.oio 


8.39 < log(MBH/MQ) < 9.0 0.476 


0 444 +011’^ 
-0.106 


n 1 -I y +0.083 
-0.036 
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^Nuo(^) 


Fig. 1.— Distribution of AGN nuclear magnitude in the /2band from K15. The red solid line 
shows the nuclear magnitude distribution of AGNs for which the host galaxy has B/T > 0.4. 
The black dashed line shows the nuclear magnitude distribution of AGNs for which the host 
galaxy has B/T < 0.4. The black dotted line represents a cut at M^uc = —22.44 mag. 
This hgure shows that the host galaxies of AGNs with Mnuc(-R) < —22.44 mag are mostly 
bulge-dominated galaxies. 
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Fig. 2.— Redshift versus absolute magnitude in the R band of the host + nuclear component 
(small circles) and the nuclear component only (Mnuc(-R); stars) for 39 AGNs. See section 
12.31 for how we derived Mnuc- 
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Fig. 3.— Redshift versus observed surface-brightness limits of 39 AGNs. The black dashed 
line is the median limit, 27.2 mag arcsec“^ at la. The red stars, the black circles, and the 
blue squares are the surface-brightness limits in the R or r, V, and i filters, respectively. 
The hlled symbols denote AGNs that have merging features (see Section H]). The dotted line 
indicates the effective surface brightness limit of the control sample images (See Section l3l^ . 
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(o) CTS J17.17 (b) Model Image (c) Residuol Imoge 



(d) PKS 0159-062 (e) Model Image (f) Residual Image 



Fig. 4.— CFHT images of two AGNs fitted by GALFIT; (a) original, (b) model, and (c) 
residual images of GTS J17.17, and (d) original, (e) model, and (f) residual images of PKS 
0159-062. The length of the horizontal bar in each panel is 10". The background of each 
image is subtracted, and the pixel values are square-rooted to bring out the low surface- 
brightness features. Identical pixel value scales are adopted to plot images for the same 
object. 
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Fig. 5.— Redshift versus (a) BH masses, (b) the host bulge magnitudes, (c) the Eddington 
ratios, and (d) the bolometric luminosities of the base AGN sample (black circles) and AGNs 
with a merging feature as classihed in Section H] (red stars). 
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Fig. 6.— Examples of the Stripe 82 early-type galaxies with merging features classihed as 
tidal tail (T), shell (S), interaction (I) types, respectively. The length of the horizontal bar 
in each panel is 10". 
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Fig. 7.— E xamples of gal axies with ambiguous classification. They are classihed as merging 
systems bv iKavirail (120101 1. However, we do not identify the merging feature, and it is likely 
that the merging feature, even if one assumes that it exists, will not be identihed as such by 
us in our AGN sample. The length of the horizontal bar in each panel is 10". 
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Fig. 8.— (a) Bulge magnitudes in the R band of the control sample early-type galaxies 
(small black circles), the control sample early-type galaxies with a merging feature (small 
red circles), AGN hosts (black stars), and AGN hosts with a merging feature (red stars) as 
a function of redshift, (b) Black hole masses of the control sample early-type galaxies (black 
circles), the control sample early-type galaxies with a merging feature (red circles), AGN 
hosts without (black stars) and with (red stars) merging feature as a function of redshift. 
The blue rectangles denote AGNs from the DuPont sample. 
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Fig. 9.— Black hole mass distribution of the control sample early-type galaxies and AGNs. 
The unhlled black histogram and the dashed line are the black hole mass distribution and 
its median of the control sample early-type galaxies. The black hlled histogram 

and the dotted line are the black hole mass distribution and its median (10®'^® Mq) of the 
AGNs. 
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Fig. 10.— Merging fraction of the control sample early-type galaxies as a function of Mbh 
and Mbui- The merging fraction increases as Mbh and Mbui increases. The horizontal bars 
show the Mbh and Mbui ranges covered by early-type galaxies in two bins divided at Mq 
and -21.81 mag. 
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Fig. 11.— Examples of simulated AGN images. The seeing is 1". The redshifts and the 
magnitudes of the nuclear point sources are varied as indicated (see text). The length of the 
horizontal bar in each panel is 10". 
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Fig. 12.— Examples of simulated AGN images. The seeing is 2". The redshifts and the 
magnitudes of the nuclear point sources are varied as indicated (see text). The length of the 
horizontal bar in each panel is 10". 
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Fig. 13.— Examples of simulated AGN images. The seeing is 1". The redshifts and the 
magnitudes of the nuclear point sources are varied as indicated (see text). The length of the 
horizontal bar in each panel is 10". 
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Fig. 14.— Examples of simulated AGN images. The seeing is 2". The redshifts and the 
magnitudes of the nuclear point sources are varied as indicated (see text). The length of the 
horizontal bar in each panel is 10". 
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Redshift 



Redshift 

Fig. 15.— Change in the merging fraction when simnlated AGNs are redshifted from z = 
0.041 to a higher redshift bnt withont adding a nnclear component (the npper panel) and 
when the simnlated AGNs are redshifted and a nnclear component is added (the bottom 
panel). The magnitnde of the nnclear component is assnmed to be mNuc = '^host — 1-0 mag. 
Gases for two types of merging featnres (black cross for T type and red hlled circles for S 
type) are indicated. The dashed and the dotted lines in the bottom hgnre denote the resnlt 
for the T type and the S type in the npper hgnre, respectively. 
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Fig. 16.— Images of 39 AGNs in our base sample. T=tidal-tail type, S=shell type, 
I=Interaction type. The length of horizontal bar in each panel is 10". 
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Fig. 16.— Continued 
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Fig. 16.— Continued 
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Fig. 16.— Continued 
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Fig. 16.— Continued 
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Fig. 17.— Merging fraction of AGNs and the simulated AGNs as a function of redshift. The 
left panel is the merging fraction of objects with log(MBH/M 0 ) < 8.39, and the right panel is 
the merging fraction of objects with log(MBH/MQ) > 8.39. The panels for each Mbh are also 
divided into three, each corresponding to a different nuclear magnitude. The magnitudes of 
the added PSF components for each panel are (a) top: mgai, (b) middle: rugai — 1.0, and 
(c) bottom: rugai — 1.5. The black asterisks are the merging fraction of early-type galaxies 
without any nuclear source (PSF) included, and the black circles are the numbers for the 
cases with a nuclear component included. The seeing of simulated images is assumed to be 
here. 
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Fig. 18.— Merging fraction of AGNs and simulated AGNs (control sample) versus the black 
hole mass of the host galaxy. The merging fraction of simulated AGNs has been corrected 
for the difference in the Mbh distribution between the AGN and the control samples. The 
red stars are for AGNs in the base sample, and the hlled red stars are for AGNs in the best 
sample (Mnuc(-R) < —22.44). The hlled black circles are for simulated AGNs. 
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Fig. 19.— Merging fraction of AGNs and simulated AGNs (control sample) versus the bulge 
magnitude of the host galaxy. The merging fraction of simulated AGNs has been corrected 
for the difference in the Mr distribution between the AGN and the control samples. The red 
stars are for AGNs in the base sample, the hlled red stars are for AGNs in the best sample 
(Mnuc(-R) < —22.44), and the hlled black circles are for simulated AGNs. 
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